Physical Consequences of Complex Dimensions of Fractals 
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It has been realized that fractals may be characterized by complex dimensions, arising from 
complex poles of the corresponding zeta function, and we show here that these lead to oscillatory 
behavior in various physical quantities. We identify the physical origin of these complex poles as the 
exponentially large degeneracy of the iterated eigenvalues of the Laplacian, and discuss applications 
in quantum mesoscopic systems such as oscillations in the fluctuation Y?{E) of the number of levels, 
as a correction to results obtained in Random Matrix Theory. We present explicit expressions for 
these oscillations for families of diamond fractals, also studied as hierarchical lattices. 
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Fractals, such as the well-known Sierpinski gasket, 
have been thoroughly studied in physics and in math- 
ematics. In addition to their own intriguing properties, 
they provide a useful testing ground to investigate prop- 
erties of disordered classical or quantum systems [1], ad- 
dressing such fundamental physical issues as Anderson 
localization, the renormalization group, and phase transi- 
tions [2] . In addition to condensed matter and statistical 
physics, fractals have been considered in other contexts 
such as gravitational systems [3, 4], and in quantum field 
theory [5] . Despite the large amount of work dedicated to 
the study of the spectra of deterministic fractals, explicit 
expressions for spectral functions such as heat kernels or 
spectral zeta functions, from which many physical quan- 
tities can be derived, have remained elusive. It is well 
known that the heat kernel Z{t) and zeta function C,{s) 
play central roles in various fields of physics: from meso- 
scopic physics [6], to black holes [7], to quantum field 
theory on curved spaces such as de Sitter and anti De 
Sitter spaces [8] , to the physics of the Casimir effect [9] . 
This is largely due to their relation to the notion of the 
partition function in statistical physics [10], and to the 
ubiquity of Schwinger's proper-time formalism [11]. 

An important step was to identify the leading contribu- 
tion to Wcyl's small time expansion of Z{t), showing that 
it is determined by the fractal's spectral dimension ds 
[12], rather than by its fractal (Hausdorff) dimension dh, 
as initially conjectured. The fact that fractals are charac- 
terized by a set of more than one dimension, as opposed 
to standard Euclidean spaces, illustrates the richness and 
peculiarity of self-similar structures. Spectral properties 
of deterministic fractals have recently been considered 
anew in mathematics, and the notion of complex valued 
fractal dimensions has been introduced [12, 13], leading 
to new results for the zeta function [14-16]. In this Let- 
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ter we use and extend these results to study the resulting 
(log-periodic) oscillations in the heat kernel and related 
physical quantities. We illustrate these ideas with a spe- 
cial class of fractals known as diamond fractals. These di- 
amond fractals permit simple explicit formulas, yet they 
exhibit properties representative of a wider class, includ- 
ing the Sierpinksi gasket (we discuss this general class in 
the conclusions). The diamond fractals also allow us to 
vary the spectral dimension, in particular to values less 
than, greater than, or equal to the critical dimension 2. 

Log-periodic oscillations have a long history in physics: 
in the theory of phase transitions [17], the renormaliza- 
tion group [18], Levy flights and fractals [19], and gen- 
erally in systems with a discrete scaling property [20]. 
Diamond fractals have also been studied in the physics 
of hierarchical lattices [21]. 

Our main result is the identification and characteri- 
zation of a new oscillating behavior of Z{t) at small t, 
which has implications for various physical quantities. 
Such oscillations do not exist for smooth manifolds, or 
even for quantum graphs. We apply these considerations 
to the concrete case of quantum mesoscopic systems [6], 
and show that the oscillating behavior can be directly ob- 
served in spectral quantities such as the fluctuations of 
the number of energy levels and the Wigner time delay. 
We also relate the electric conductance g, the associated 
weak localization corrections lS.g, and universal conduc- 
tance fluctuations 5g^ to the fractal zeta function. 

We first recall some basic definitions and facts about 
deterministic fractals. As opposed to Euclidean spaces 
characterized by translation symmetry, self-similar (frac- 
tal) structures possess a dilatation symmetry of their 
physical properties, each characterized by a specific frac- 
tal dimension. To illustrate them, we consider through- 
out this letter the family of diamond fractals (see Fig. 
1), but keeping in mind that our results apply to a much 
broader class of fractals, including the Sierpinski gasket. 
At each step n of the iteration, we characterize a fractal 
by its total length L„, the number of sites A^„, and the 
diffusion time T„. Scaling of these dimensionless quan- 
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titles allows to define the corresponding Hausdorff dh, 
spectral ds, and walk d^ dimensions according to 



dh = 



In Lr, ' 



dw = 



In Lr, 



ds = 2 



In Tn 



(1) 



where the limit n — > oo is understood. These three di- 
mensions are thus related by d^ — 2dh/dw 




For example, from the spectral determinant, we deduce 
the density of states: p{E) = —■^ limj^o+ Im^ In 5(7), 
with 7 = —E + ie. This can also be written [23] in 
terms of the on-shell S-matrix S{E) by the Birman-Krein 
formula p{E) = In det iS(— i?), also defining the 

Wigner time delay: t{E) = — In det S{-E). 
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FIG. 2: Fractal dimensions and size scaling factor for diamond 
fractals and for the Sierpinski gasket. For -De, 2, the spectral 
dimension is ds ~ 2.58, and for De.a, ds ~ f.63. 




FIG. 1: First 2 iterations of the diamond fractals -D4,2, -De, 2 
and -De, 3. Their respective branching factors (defined in the 
text) are -B = f,2, 1. 

To obtain the heat kernel of a fractal, let us recall 
the corresponding expression for an Euclidean system 
of space dimension d. We consider the diffusion equa- 
tion —Atl'k{r) = -Efc-0fe(r), where the diffusion coefficient 
is set to unity, without yet specifying boundary condi- 
tions. The probability P(r,r',t) to diffuse, in time t, 
from an initial point r to a final point r', is given by 
the Green's function defined in an arbitrary volume fJ: 
P{vy,t) = e(OEfeE,ti^fe,,(r)^fe,j(r')e-^^*. Here 
gk is a degeneracy factor generally different from unity 
{e.g. on a sphere [22]), except for one dim. diffusion on a 
finite interval. The heat kernel Z{t) is defined for t > 0: 



Z{t) = f P(r,r,<)dr = Vfffee- 
Jn 



■E,.,t 



(2) 



The spectral zeta function is defined by a Mellin- 
Laplace transform of the heat kernel 



C(s,7) = 



1 
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gk 



{Ek + 7)^^ 



(3) 



Many quantities are derived directly from the spectral 
zeta function. E.g., the spectral determinant 5'(7) is [7] 



S'(7) = det(-A + 7) = exp 



ds 



C(s,7)l 



(4) 



which follows directly from the analytic continuation of 
C(s, 7) in the complex s plane as a meromorphic function 
analytic at s = 0, and the identity [^A~*] 



s=0 



In A. 



To generalize (2) to a fractal, we consider the proba- 
bility P(r, t) to diffuse over a distance r in a time t (with 
obvious notations). Scaling properties of diffusion are 
expressed using the definition (1) of the walk dimension 
dw through the scaling transformation, P{Xr, X^^'t) — 
P{r, t), for any scaling factor A of the length, so that the 
probability is of the form P{r,t) = f{r'^'"/t), where / is 
some unknown function. In addition, the normalization 
condition, J d'^'^rP{r, t) = 1, and the change u = r/i^/''™, 
lead to the general scaling form 



P(r, t) 



tdh/d. 



(5) 



This implies that diffusion on a fractal is anomalous in 
the sense that the usual Euclidean relation (r^(t)) cx t, 
for long enough times, is now replaced by (r^ (t)) cx t^/''™ : 
hence the name "anomalous random walk dimension" 
for dw. Then, relations (1) imply the well known re- 
sult P(0,t) oc for the leading term of the return 
probability which is driven by the spectral dimension ds , 
rather than by the Hausdorff dimension dh. Generalizing 
(2), the heat kernel of a diamond fractal can be obtained 
by noticing that the spectrum of diamond fractals is the 
union of two sets of eigenvalues. One set is composed of 
the non degenerate eigenvalues Tr^fc^, (for k = 1,2,.. .). 
This corresponds to the spectrum of the diffusion equa- 
tion defined on a finite one-dimensional interval of unit 
length, with Dirichlet boundary conditions. The sec- 
ond ensemble contains iterated eigenvalues, Tr^fc^L^™ , ob- 
tained by rescaling dimensionless length Ln and time Tn 
at each iteration n according to L^^ = Tn, given in (1). 
To proceed further, we use the explicit scaling of the 
length Ln = I" upon iteration (see Table). These iterated 
eigenvalues have an exponentially large degeneracy given, 
at each step, by SLf,'- = B {l'^")", where B = (l'^"-^ - 1) 
is the branching factor of the fractal (see Fig.l), and the 
integer l'^'^ is the number of links into which a given link 
is divided. The exponential growth of the degeneracy 
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plays a crucial role in our analysis. By contrast, on an 
-/V-dimensional sphere the degeneracy grows as a poly- 
nomial, of order A'' — 1 [22] . Finally, the diamond heat 
kernel Zjj^t) is the sum of contributions of the two sets 
of eigenvalues: 



Jeadin 



oo oo 



Zoit) = e-^^ * + bY, Lt E e-'^^ . (6) 

k=l ri=0 k=l 

The associated zcta function Cr>(s), from (3) at 7 = 0, is 



Cd(s) 



\ n=0 / 



(7) 



where ^ji{2s) is the Riemann zeta function. Note that 
a very similar structure arises for the Sicrpinski gasket 
[14], with the Riemann zeta function factor replaced by 
another zeta function. Cd{s) has complex poles given by 
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(8) 



where m is an integer. The origin of these complex poles 
is clearly the exponential degeneracy factors. The com- 
plex poles have been identified with complex dimensions 
for fractals [13, 14]. 

By an inverse Mellin transform, wc can write the heat 
kernel as Zoit) = /;+^ CD(s)r(s) ^-^ Then the 
leading small time behavior comes from the pole of Cd(s) 
at s = sq = rfs/2, giving the anticipated time decreasing 
function i"'''/^. The pole of Cd(s) at s = 1/2 (coming 
from the C-r(2s) factor) has zero residue for all diamonds, 
and so does not contribute to the short time behavior 
of Zd (t) . (Remarkably, this vanishing of the residue at 
s = l/dw also applies to the analogous zeta function on 
the Sierpinski gasket [14]). The pole of T{s) at s = 
gives a constant contribution, C-d(O), to Zjy{t). But the 
really surprising new behavior comes from the complex 
poles in (8), leading to the oscillatory behavior: 



Zoit) 



dh-l 



- 1 1 



Oo + 2Rc \ ait 
+ Cl3(0) + . 



-2j:ir/(d„ In;) 



(9) 



where we have defined a™ = r(sm)C-R(2sm)/7''^*''" ■ The 
leading term oc t~'^'l'^ is therefore multiplied by a peri- 
odic function of the form a\r cos(lnf;''i') + an sin(lnt*i'), 
where a\r,i are respectively the real and imaginary parts 
of ai, and su = 27r/lnZ'*™. The oscillations of Zo{t) 
are represented in Fig. 3, and we note that the higher 
complex poles give much smaller contributions. Similar 
behavior has been found numerically for the Sierpinksi 
gasket [24]; from our work, we further find explicit ex- 
pressions for the coefficients, also in the Sierpinksi case. 
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FIG. 3: Heat kernel Zoit) at small time, normalized by the 
leading non-oscillating term, for the fractal diamond 7?4,2. 
The solid [blue] curve is exact; the dashed [red] curve is the 
approximate expression (9). At very small t, these curves are 
indistinguishable, as shown in the inset plot. The relative 
amplitude of the oscillations remains constant as t — > 0. 



In principle, all spectral properties can be derived from 
the heat kernel (6), or from the associated zeta function 
Cd{s) in (7), even though those are not directly accessi- 
ble physical quantities. For example, the constant term 
Cd{0) in (9) leads to a topological term (d{0)S{E) in the 
density of states. More interestingly, the oscillations of 
Zo{t) lead to oscillatory behavior in physical quantities. 

We give an explicit example of one such quantity, in 
quantum mesoscopic systems. The fluctuation T,'^{E) of 
the number of levels within an energy interval of width E 

is defined by the variance, T,'^{E) = 7P(E) - N{Ef, of 
the integrated density of states (the counting function). 
In the diffusion approximation, one can express S^(i?) 
directly in terms of the heat kernel through [6] 



T.'^E) 



dt 



Zoit) sin^ 



Et 
T 



(10) 



Inserting (9) for Zo(t), we obtain. 



Y.^{E) 



idH-i _ I 
2tt In 



S'*-'/^ f6o + 2Rc ffei^^.vm/'' 



idH- 



27r2(?'ih _ 1) 



{I -l)lii4:E -I 



Idu-i _ I 
l'^'^ - 1 



hif- (11) 



where = Cj?,(2sm)/(srn7r^'''" sin(7rsm)). The leading 
term oc E"^"^^ is now multiplied by a periodic function 
of the form 61^ cos(ln£^'*") -|- 6ii sin(lni?*"), where bir,i 
are respectively the real and imaginary parts of bi . This 
oscillating behavior of Tj'^{E) is represented in Fig. 4. It 
is remarkable that the behavior of ^^(i?) differs at low 
energy from the expected ergodic regime independent of 
fractal dimensions which is well described by Random 
Matrix Theory, and also at large energy from free diffu- 
sion oc iJ'^s/^ on the fractal. 
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FIG. 4: Semi-log plot of the fluctuation T?{E), normalized by 
the leading non-oscillating part, for the diamond fractal 1)4,2. 
The solid [blue] curve is exact; the dashed [red] curve is the 
approximation (11). At very large E, these curves are indis- 
tinguishable, as shown in the insert. The relative amplitude 
of oscillations remains constant as £ ^ 00. 



results apply. We have chosen to illustrate our results 
using the class of diamond fractals because this class 
permits simple and explicit formulas. But the general 
observations about complex dimensions of fractals, and 
oscillations in the heat kernel trace and associated physi- 
cal quantities, generalize to the class of fractals known as 
finitely ramified self-similar fractals with full symmetry 
group (i.e. the symmetry group has doubly transitive ac- 
tion on the boundary). A complete description of these 
fractals can be found in [16, 29, 30]. The best known ex- 
amples are the Sierpinski gasket ([24]), the Level-3 Sicr- 
pinski gasket and the Vicsek set. Note that in these and 
the majority of other examples the walk dimension is not 
2 [as it is for all diamond fractals], and so dg is not nec- 
essarily equal to d^,. 
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Transport quantities such as the dimensionless conduc- 
tance (expressed in units of e^/h), are also interesting 
on a fractal. For instance, the so-called weak localiza- 
tion correction to the conductance, and conductance 
fluctuations described by the variance 5g^, a universal 
quantity independent of the system size, take a general 
and remarkable form [6] expressed in terms of the zeta 
fimction only, namely Ap = — 2C_d(s = 1) [N.B. the pole 
must be treated properly in the special case ds = 2], and 
5g^ = \2C,d{s = 2). For a single channel setup, it is also 
possible to relate the Fano factor [25] which characterizes 
shot noise to A5 [26] and thus to the zeta function, and 
we obtain immediately F = — Ag = 1/3 and 5g^ — 2/15 
for diffusion on a finite one-dimensional interval. On a 
fractal, and using (7), these quantities now depend on the 
fractal dimensions dh and ds, and therefore conductance 
experiments could be used to determine them. 

To summarize, we have considered spectral properties 
of deterministic fractals such as the heat kernel and the 
spectral zeta function. Using the class of diamond frac- 
tals, we have derived simple and explicit formulas which 
illustrate a new and general oscillatory behavior of the 
heat kernel, and relate it to complex poles, also identified 
as complex fractal dimensions, resulting from the expo- 
nentially large degeneracy of the iterated eigenvalues of 
the Laplacian. These oscillations which show up in a 
variety of interesting physical quantities, characterize a 
fractal. Our results may be useful to study properties 
of more general quantum graphs [6] where degeneracies 
must properly be taken into account, and to investigate 
magnetic [27] and topological properties [28] of fractals 
when submitted to external fields such as a Aharonov- 
Bohm fluxes. They may also have interesting implica- 
tions for gravitational and quantum field theoretic appli- 
cations. 

We end by stating the class of fractals to which our 
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